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Abstract: In the work, transmission of nerve impulses along nerve fibers is simulated. Research is being conducted in 
the framework of the electrical theory of propagation of the action potential. The soliton approach is used. The numerical 
experiment on the evolution of the nerve impulse and the laws of motion is conducted. It is shown that the stable form of 
nerve impulses solitons is realized under different initial conditions. The threshold character of occurrence of nerve pulse is 
simulated. It is shown that the number of pulses produced changes depending on the degree of nonlinearity: with strong 
initial excitation of large amplitude soliton is unstable; it breaks up into a multitude number of solitons with small ampli- 
tudes. In this case, the greater the nonlinearity parameter, the greater the number of births of solitons. Unusual solitonlike 
regimes of interaction of nonlinear pulses excitations are illustrated; in certain anomalies, colliding nerve impulses are re- 
flected instead of their usual quenching. The possibility of the decay of the nerve impulse at the bifurcation of the nerve 
fibers or the appearance of heterogeneity of the passing of a nerve impulse (the presence of dissipation in the environment) 
is investigated. Physical modeling allows obtaining and studying the analytical results and elucidating the physical prin- 
ciples of biological processes. 
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1. Introduction 

To coordinate work of the different systems in the hu- 
man body requires harmonious work of its most complex 
system - the nervous system. Its basic function is the prop- 
agation of an electric pulse through the nerve fiber. Physio- 
logical studies have shown that propagation of a nervous 
pulse through neural nets has auto wave character. T. 
Heimburg and A.D. Jackson proposed a new neural theory 
called Soliton theory. In this theory, the nerve conduction 
is proposed as a density wave. [1] B.Lautrup, R.Appali 
and others [2] examined a class of solitons, obtained from 
a generalization of the Boussinesq equation. Soliton model 
is based on the propagation of a localized density wave in 
the axon membrane. This allows simulating such processes 
using Soliton Theory, which describes the form and veloci- 
ty of propagation of the nervous pulse, stability of its cha- 
racteristic parameters in time, and process locality with 
perfect precision. [3] 

According to electrical theory, the process of transmis- 
sion of the pulse along the nerve fiber is a process of prop- 
agating the action potential: the nerve fiber extends the 
electrochemical reaction, which generates a running im- 
pulse of voltage. The nerve pulse or action potential propa- 



gates through the nerve fiber without decay and with con- 
stant velocity. An impulse can be formed and spread be- 
cause the nerve fiber contains a nonlinear element, which 
suppresses small deviations from normal state and streng- 
thens big ones. In this case, the nonlinear dependence of 
the membrane permeability (through which the nerve im- 
pulse is spread) from the momentum and the diffusion of 
ions across the membrane are balanced. In the paper [4] is 
shown that nonlinearity can actually balance diffusion and, 
as a result, a running solitary wave with constant speed and 
shape can occur, i.e. a soliton. Thus a nerve impulse is an 
extension of an action potential. Physiological studies have 
shown that propagation of a nervous pulse through neural 
nets has autowave character. This also allows simulating 
such processes using Soliton Theory which describes the 
form and velocity of propagation of the nervous pulse, sta- 
bility of its characteristic parameters in time, and process 
locality with very good precision. [5] 

2. Theory and Numerical Experiment 

Since the nerve impulse is resistant, it extends without 
distortions and significant loss of energy; its movements 
can be modeled as a Korteweg de Vries soliton. In this case 
it is described by the nonlinear differential equation [6]: 
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u(x, 0) = u 0 exp 




(3) 



where the second term characterizes a weak (cubic) dis- 
persion, and the third term describes nonlinearity. The 
solution of this equation has the following form, 



U2( X ~ Vt \ 

ch 

where u is the action potential , u max is the amplitude of 
the action potential (in our case of the nerve pulse), v is the 
velocity 



(v Umax/ 3) (2 Cl) 

, and A is determined from the ratio 

A= V 12/3 Aw, (2 b). 

From (2 a) it follows, that the propagation velocity of a 
pulse increases with increase in the amplitude, whereas 
from (2b) it follows that the soliton width decreases with 
the increase in the amplitude. Therefore, the solitons line 
up in the order of decreasing amplitudes and growing 
widths “Fig. 1”. 





Figure 1 . Soliton distributions at different times. (Numerical experiment). 

We performed a numerical experiment of formation, 
evolution and movement mechanism of solitons. The stable 
form of solitons can be obtained with different initial con- 
ditions. 

We chose the initial perturbation in two forms: 



(where / is the linear size of initial perturbation) and a 
hyperbolic cosine. However, regardless of the character of 
initial perturbation, over time the pulses transform into a 
soliton. Obviously, the initial conditions are connected 
with the signal received from the brain through conducting 
ways and depend on the state of the organs or the mental 
state of the person. ( These data are consistent with [6] ) 
English physiologists A. Hodgkin and A. Huxley have 
developed a theory (based on numerous experiments) 
where they concluded that the shape and the speed of an 
impulse don’t depend on the stimulus of the nerve. If the 
stimulus is high, several consecutive impulses are produced, 
while if it is low, no impulse will be produced through the 
nerve; the minimal power of stimulus to produce an 
impulse is called threshold. Each impulse transfers one unit 
of information and the receiving devices count the number 
of elementary particles of information received and the 
time it took [7]. All this reminds of the movement of 
solitons. In paper [8] Hodgkin-Huxley and Soliton theories 
are described and a theoretical comparison has been carried 
out throughout the analysis of the theories and models. 

In the Korteweg de Vries equation the degree of nonli- 
nearity is characterized by a dimensionless parameter 
a = l(u 0 / fl) 1 / 2 . The numerical experiment showed that at 
a < (J m iri the solitons are absent, and instead sine wave 
oscillations arise “Fig. 2”. 
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Figure 2. Sine wave oscillations. o < o mil 
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From “Fig 2a” it is visible that solitons are absent, but to 
the left of an initial impulse there are oscillations. This case 
is realized at /3 > 0. At /3 < 0 (“Fig 2b”) oscillations 
arise on the right. The received result follows from Korte- 
weg de Vries equation. 

We investigate character of the solution of the equation 
(1) depending on values a. For a = 5.9 we studied evolu- 
tion of an initial impulse for times t x < t 2 < t 3 . Thus the 
initial impulse moves to the right, the first solitons ampli- 
tude grows over time, at t 2 the second soliton is born 
whose amplitude also grows over time; at t 3 we have 2 
solitons, and the amplitude of the first one is higher than 
the amplitude of the second one. At a = 9 and the same 
initial condition the decision for the same values of time is 
received. It can be seen from the figure that 2 separate soli- 
tons on the right and the third, a new one, on the left are 
visible. Thus, with increase of a at an identical initial im- 
pulse the number of born solitons increases “Fig. 3”. 

It is known that in vivo (living organism) nerve impulses 
moving along the axon toward each other collide and anni- 
hilate each other, which provides one-way transmission of 
signals from receptors in the central nervous system. How- 
ever, the works of Aslamanidi O.V., Momev O.A. [9,10] 
are based on numerical experiments with Hodgkin-Huxley 
equations predict a new effect: for some values of the pa- 
rameters colliding nerve impulses may reflect. It is shown 
that in the soliton mode, the propagating impulse is a 
“doublet” consisting of a high-amplitude impulse-leader 
and an accompanying low-amplitude wave. The first colli- 
sion of the leaders causes annihilation, and the second col- 
lision of the low-amplitude waves leads to summation. As 
a result, at the point of collision the potential reaches the 
above-threshold value, causing regeneration of “doublets” 
that become scattered in different directions. 
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Figure 4. Collision of two solitons. The figures show successive moments 
of time. 



Figure 3. The different number of originating solitons for different values 
of nonlinear parameter. 

We have modeled the collision of two solitons occurring 
in a similar mode, on the basis of Korteweg de Vries equa- 
tion “Fig. 4”. 



u 




The analysis of such regimes of interaction of nonlinear 
pulse can be useful in deciphering the mechanisms of car- 
diac arrhythmias. It is known that the Korteweg de Vries 
soliton gradually slows down under friction, reducing 
height and becoming blurred. The lifetime is inversely pro- 
portional to the force of friction. Analogically, when the 
cross section is changed or nerve fibers are bifurcated, the 
passage of nerve impulses may become difficult or even 
completely blocked. In the equation, this can be considered 
as the presence of friction. [6] We have also investigated 
the solitons’ evolution in the presence of friction. When 
taking friction into account, the Korteweg de Vries equa- 
tion becomes the form 



du du d 2 u d 3 u 
at +U 5 t ~ U 'dt I + ^dx^ = ° 



ti 



u 




t 2 



u 




where ft is the coefficient of friction. Suppose that at 
t = 0 we have an initial perturbation in the form 

u(x, 0) = u 0 ch _2 (x/0- 

At (i = 2 the numerical experiment was made for t= tl< 
t2< t3< t4< t5. “Fig. 5a”. Next, the same experiment was 
performed for (i = 10 and for the same values of time. As 
expected, attenuation occurs faster with with the increasing 
coefficient of friction. “Fig.5b”. 
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Figure 5. Solitons evolution at different coefficients of friction ju for the 
same values of time. 

3. Conclusion 

Thus, different regimes of propagation of nerve pulses 
can be realized under different conditions. Solitonlike solu- 
tions describe a wide range of phenomena that accompany 
the transmission of nerve pulses, heart rate, self-oscillation. 
The soliton models of processes make it possible to obtain 
adequate description and identification of the state of such 
systems as the conducting muscular and circulatory sys- 
tems of the heart. [11] By changing parameters of the soli- 
ton model, it is possible to model the pathologies of nerv- 
ous conducting ways of the vascular and muscular system 
of the heart. The soliton approach opens up large possibili- 
ties in the research into the complex processes in living 
organisms. Possible pathology in the vascular system 
changes spectral picture of the soliton waves and allows 
diagnosing a variety of diseases of nerve pathways and the 
circulatory system. Mathematical modeling provides the 
analytical results and the way to identify the physical prin- 
ciples of biological processes. 
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